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Abstract 



Phantom dark energy is a proposal that explains the current observations that mildly favor the equation 

of state of dark energy oj'*'^ crossing —1 at 68% confidence level. However, phantom fields are generally 

ruled out by ultraviolet quantum instabilities. To overcome this discrepancy, in this paper we propose a 

mechanism to show that how the presence of bulk viscosity in the cosmic fiuid can temporarily drive the fluid 

into the phantom region (to < —1). As time is going on, phantom decays and ultimately cj'*'^ approaches to 

\^ ' —1. Then we show these quintessence and phantom descriptions of non-viscous and viscous dark energy and 

(-H . reconstruct the potential of these two scalar fields. Also a diagnostic for these models are performed by using 

Oj' the statefinder pairs {s,r}. All results are obtained in an anisotropic space-time which is a generalization of 

I ' FLRW universe. 

o: 
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1 Introduction 



^T^ [ It is a very well known fact that our universe is experiencing an accelerating expansion at the present time (Perl- 
Si ■ mutter et al. 1999; Riess et al. 1998, 2001; Tonry et al. 2003; Tegmark et al. 2004). It is believed that an exotic 
psj [ form of energy with negative pressure called dark energy is responsible for the current observed accelerating 
.• ■ expansion of the universe (Tegmark et al. 2004; Bennet et al. 2003; Spergel et al. 2003; Abazajian et al. 2004). 
r— s I According to the recent observations we live in a nearly spatially flat Universe composed of approximately 4% 
^f) ■ baryonic matter, 22% dark matter and 74% dark energy. However, the observational data are far from being 
complete. It is not even known what is the current value of the dark energy effective equation of state (EoS) 
parameter w'^'^'^^ — p(''^)/p('''^) which lies close to —1: it could be equal to —1 (standard ACDM cosmology), a 
little bit upper than —1 (the quintessence dark energy) or less than —1 (phantom dark energy). One of the main 
KN ■ candidate for dark energy is cosmological constant A, which has pressure p^'^'^-' = —p'^'^''\ Although, cosmological 
i3 I constant can explain the current acceleration phase of universe, it would suffer from many serious theoretical 
" ■ ■ problems, such as the fine-tuning and the coincidence problems. Another candidate for dark energy is provided 
by introducing scalar fields. An important class of scalar fields are known as "quintessence" with — ^ > w > — 1 
(Ratra and Peebles. 1988; Wetterich 1988; Turner and white 1997; CaldweU et al. 1998; Liddle and Scherrer 
1999; Steinhardt et al. 1999) in which the scalar field mimics the perfect fluid and hence could lead to a solution 
for coincidence problem. However, quintessence scenario of dark energy is not in accurate consistent with recent 
observations as w < —1 has been favored by recent observations (Knop et al. 2003; Riess et al. 2004; Alam et al. 
2004; Hannestad and E. Mortsell 2004). To get w < — 1, a new class of scalar field models with negative kinetic 
energy, known as "phantom field" models have been suggested (Caldwell 2002). Nevertheless, in this case the 
universe shows some very strange properties (Carroll et al. 2003; Cline et al. 2004; Buniy and Hsu 2006; Buniy 
et al. 2006). For example, since the energy density of phantom field is unbounded from below, the vacuum 
becomes unstable against the production of positive energy fields hence these fields are generally ruled out by 
ultraviolet quantum instabilities (Carroll et al. 2003). Another problem is the future finite singularity called Big 
Rip (Caldwell et al. 2003) which leads to the occurrence of negative entropy (Brevik et al. 2004). Therefore, on 
the one hand observations mildly favors models with uj crossing —1 near the past and on another, models with 
w < — 1 are unstable from theoretical point of view. In this paper we suggest a simple mechanism to overcome 
this discrepancy by introducing bulk viscosity in the cosmic fluid. First, viscosity causes dark energy which is 



varying in quintessence to pass the phantom divided line (PDL) and drop it to phantom region. Next, since 
viscosity is a decreasing function of time, it will die out and w will leave phantom region and tend to —1 at late 
time. Hence the problem of future singularity (big rip) will never occur in this scenario. 

It has been shown in refs (Mclnnes 2002; Barrow 2004) that, an ideal cosmic fluid, i.e. non-viscous, give 
raise to the occurrence of a singularity of the universe in the far future called big rip. The singularity problem 
can be modified or soften via following two methods. The first is the effect of quantum corrections due to the 
conformal anomaly (Brevik and Odintsov 1999; Nojiri and Odintsov 2003, 2004) and second, is to consider the 
bulk viscosity of the cosmic fluid (for example see (Misncr 1968; Padmanabhan and Chitre 1987; Brevik and 
Ballanger 2004). The viscosity theory of relativistic fluids was first suggested by Eckart, Landau and Lifshitz 
(Eckart 1940; Landau and Lifshitz 1987). The introduction of viscosity into cosmology has been investigated 
from different view points (Gr0n 1990; Barrow 1986; Zimdahl 1996; Maartens 1996). The astrophysical observa- 
tions also indicate some evidences that cosmic media is not a perfect fluid (Jaffe et al. 2005), and the viscosity 
effect could be concerned in the evolution of the universe (Brevik and Gorbunova 2005; Brevik et al. 2005; 
Cataldo et al. 2005). It was also argued that a viscous pressure can play the role of an agent that drives the 
present acceleration of the Universe (Zimdahl et al. 2001; Balakin et al. 2003). The possibility of a viscosity 
dominated late epoch of the Universe with accelerated expansion was already mentioned by Padmanabhan and 
Chitre (Padmanabhan and Chitre 1987). Brevik and Gorbunova (2005), Oliver et al (2011), Chen et al (2011), 
Jamil and Farooq (2010), Cai et al (2010), Setare (2007a, 2007b, 2007c), Sctare et al (2007), Setare and Sari- 
dakis (2009), Setare et al (2009), Amirhashchi et al (2011 a, 2011b, 2011 c, 2013), Pradhan et al (2011a, 2011b, 
2011c), Saha et al (2012), and Sheykhi and Setare (2010) have studied viscous and non-viscous dark energy 
models in different contexts. Recently, viscous dark energy and generalized second law of thermodynamics has 
been studied by Setare and Sheykhi (2010). 

To be general, we use generalized FLRW equations by considering an anisotropic metric as the line-element of 
the universe. The reason for this choice of metric is behind the fact that because of high symmetry, FLRW models 
are infinitely improbable in the space of all possible cosmologies. The high symmetry involved in FLRW models 
requires a very high degree of fine tuning of initial conditions which is extraordinary improbable. Moreover, we 
can always ask that does the universe necessarily have the same symmetries on very large scales outside the 
particle horizon or at early times? 

The plan of our paper is as follows: In section 2 we give the metric and field equations. In section 3 we drive 
the generalized FLRW equations by solving the field equations of section 2. The general form of non- viscous 
and viscous dark energy equation of state parameter EoS are given in section 4. We suggest a correspondence 
between the non- viscous and viscous dark energy scenario and the quintessence and phantom dark energy model 
in section 5. In section 6, a statefinder diagnostic has been presented. In section 7 we apply our general results 
to a toy model in order to test the proposed mechanism. Our results are summarized in section 8. 

2 The Metric and Field Equations 

We consider the Bianchi type I space-time in the orthogonal form as 

ds^ = -dt^ + A^{t)dx^ + B^{t)dy^ + C^{t)dz'^, (1) 

where A{t), B{t) and C{t) arc functions of time only. 

The Einstein's field equations ( in gravitational units 8ttG = c = 1) read as 

R^-^Rg}=TJ"'''' + TJ''^\ (2) 

where T*'™ and T '^ are the energy momentum tensors of barotropic matter and dark energy, respectively. 
These are given by 



TJ"^* - diag[-p("),p("),p("),p(")], 

= diag[-l,a;("),a;("),a;(™V", (3) 



and 



rf^)'=diag[-p('^^),p('^^),p('^^),p('^^)] 



where p^"^^ and p^™^ are, respectively the energy density and pressure of the perfect fluid component or matter 
while w^™^ — p'^"^'^ / p{m) is its EoS parameter. Similarly, p^'^^'> and p^'*^) are, respectively the energy density and 
pressure of the DE component while w^'*'^^ = p(''<^)/p('*'^) is the corresponding EoS parameter. We assume the 
four velocity vector u* = (1, 0, 0, 0) satisfying u^Uj = —1. 

In a co-moving coordinate system (u* = (5q), Einstein's field equations © with ([3]) and ^ for B-I metric ([l} 
subsequently lead to the following system of equations: 

l + § + l§ = ~""'^"-"'''^''' (^) 

I + l + ll = -""^"^"''^''' (^) 

-A_D /lO tSLy J . 

If we consider a = [ABC)^ as the average scale factor of Bianchi type I model, then the generalized mean 
Hubble's parameter H defines as 

.. a I ( A B C\ 

a i \A B C j ^ ' 

The Bianchi identity Gf^ = leads to Tf- — 0. Therefore, the continuity equation for dark energy and 
baryonic matter can be written as 

p" + 3H{1 + u;"')p"' + /'^ + 3H{1 + a;'''=)p'*'= = 0. (10) 

3 Friedmann-Like Equations 

In this section, we derive the general solution for the Einstein's field equations ©-([HI). 

Subtracting Eq. (O from Eq. dH), Eq. ^ from Eq. (O, and Eq. ^ from Eq. ^ we obtain 

A B C [A b\ ^ 

a-b + c[a-b)-'^ (^^) 



B C A B C 
B^C^A\B^C^ 

and 

A C B ( A C^ 

'A^C^~B\~A^C^ 
First integral of Eqs. ^Q, ^^ and (US]) leads to 



and 



0, (12) 



(13) 
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(14) 

(15) 
(16) 



where fci, k2 and k^ are constants of integration. By taking integral from Eqs. ^^, ^T^i and (ITS)) we get 

^ = diexp[ki [{ABCy^dt], (17) 



and 



d^expiks HABCy^dt] 



A 
C 



(18) 
(19) 



where, di, ^2 and da are constants of integration. 

Now, we can find all metric potentials from Eqs. (|17p . ([T9l) as follow 



and 

Here 

ai = (did2)3, a2 = {d^^ds)^, 
where 



v4(<) = aio exp{bi / a dt), 
i3(i) = a2a exp{b2 / a^ dt), 

C{t) ~ a^a expibj, I a^ dt). 



(20) 
(21) 

(22) 



03 = (^2^3) ^ , 



fci + A:2 



fc3 - fci 



^2 + ^3 



010203 = 1, 6;^ + 62 + &3 = 0. 

Therefore, one can write the general form of Bianchi type I metric as 



ds^ = -dt 



2 ' a" 



^2g26i / a-\lt^^2 ^ ^2g26. / a-'dt^y2 ^ ^2^263 / a'^dt^^S 



Using eqs. (I^(I1) - (P^ in eqs. ([S])-® we can write the analogue of the Friedmann equation as 

2 



P 



Ka- 



and 



2(^)^> + 3.). 



(23) 



(24) 



(25) 



Here p = p™ + p'^'^, p = p™ +p'^'^ and A' == 6162 + ^163 + i'2&3- Note that iiT denotes the deviation from isotropy 
e.g. K = Q represents flat FLRW universe. Thus, when the universe is sufficiently large, almost at the present 
time, the space-time ([T]) behaves like a flat FLRW universe. 

4 Dark Energy Equation of State 

In this section we obtain the general form of the equation of state (EoS) for the viscous and non viscous dark 
energy (DE) uj'^'^ in Bianchi type I space-time when there is no interaction between dark energy and Cold Dark 
Matter(CDM) with 0;^ = 0. In this case the conservation equation (|10p for dark and barotropic fluids can be 
written separately as 

/•= + 3iJ(l + a;'*'^)/^ = 0, (26) 



and 

Eq.dlZl) leads to 



p™ + 3iJp" = 0. 



p" = p^a-\ 



(27) 
(28) 



Using eqs. (|24l) . (1^51) in eqs. (O, (|5]) we obtain the energy density and pressure of dark fluid as 

pde ^ 3^2 _ 3^^-6 _ ^m^-3 (29) 

and 

pde ^ _2- - i/2 _ ia-6^ (30) 

a 

respectively. Therefore, the equation of state parameter (EoS) of DE in it's general form is given by 



^ ^ 3 + 3La-6i/-2 - 3f7™a-3 ' 






where g = — ^^^ is the deceleration parameter, fi™ is the current value of matter density and L = 62 + ^i + ^2^3 
is a positive constant (Note that K + L = Q). 

From eq. ((3T|) we see that at present time (i.e L — 0, q = —0.55, 17™ = 0.24, a = 1), approximately, 
Wp? = —0.92. At late time, EoS parameter is given by 

u^tf - ^> (32) 

here subscript 'p/' refers to "perfect fluid" . 

According to the observations deceleration parameter is restricted as — 1 < g < 0. Therefore, from eq. ([SS]) 
we observe that at the best approximation the minimum value of ojfj? is —1 i.e EoS of non- viscous DE can not 
cross phantom divided line (PDL). In another word, non-viscous dark energy can be described by quintessence 
{uj'^^ > —1) rather than phantom {uj'^'^ < —1) field. In other hand models with w'^^ crossing —1 near the past 
have been mildly favored by the analysis on the nature of dark energy from recent observations (for example see 
(Astier et al. 2006)). SNe la alone favors a ut larger than —1 in the recent past and less than —1 today, regardless 
of wether using the prior of a flat universe (Alam et al. 2004; Astier et al. 2006) or not (Dicus and Repko 2004). 
In what follows we show that the possibility of crossing PDL will be available in a universe dominated by viscous 
DE instead of perfect DE. 

In Eckart's theory (Eckart 1940) a viscous dark energy EoS is specified by 

pt}=p';,}+U. (33) 

Here 11 = —S,{p'^'^)u\ is the viscous pressure and H = -^ is the Hubble's parameter and subscript 'f /' referees 
to "viscous fluid" . On thermodynamical grounds, in conventional physics f has to be positive. This is a con- 
sequence of the positive sign of the entropy change in an irreversible process (Nojiri and Odintsov 2003). In 
general, ^{p''"'^) — (.nip'^'^Y , where ^0 > and r are constant parameters. Note that, here we have to assume 
T > since for negative r this form of bulk viscosity does not allow our models to cross PDL. A power-law 
expansion for the scale factor can be achieved for t = 5 (Barrow 1987, 1988). It has been shown by Goliath and 
Ellis (1999) that some Bianchi models isotropise due to inflation. 



Substituting eq. p3p in eq. pip by considering the above description we obtain the EoS parameter of viscous 
DE as 

, .de _^P'' . li _ 2q-l- La-^H-^ gi-^" 

vf pde ^ pde 3 + iLa-^H--^ - 3f7S^a-3 ^^ (3f7'^«)" ' ^ ' 

, de 

where ft — jjjj and a = I — t. 

From eq. (|34p we observe that the EoS of viscous DE at present time (i.e L = 0, q = —0.55, Hq = 70, fJ™ = 

0.24, JIq"^ — 0.76, a = 1), approximately is 



213Co 
(12501.68)"' 



<^-o.92- ,^,:::^i,^ , (35) 



which clearly cross the PDL for appropriate values of a and ^o- As mentioned before, phantom flelds are 
generally plagued by ultraviolet quantum instabilities. Naively, any phantom model with uj'^'^ < —1 should 



decay to uj'^'^ = — 1 at late time. As mentioned in (Carroll et al. 2003), this ensures that there is no future 
singularity (Big Rip); rather, the universe eventually settles into a de Sitter phase. Here we highlight since 
£,{p'^'^) = S.oip'^'^y , and p'^'^ is a decreasing function of time in an expanding universe we conclude that the bulk 
viscosity dies out as time goes on and viscous phantom DE is an unstable state (as expected) and EoS of DE 
tends to —1 at late time (de-Sitter Universe). 

5 Correspondence Between Dark Energy And Scalar Fields 

It is believed that the current accelerated expansion is driven by a dynamical scalar field 4> with potential V{(j)). 

These models introduce a scalar field that is minimally coupled to gravity. As it is shown in previous section, 

one can generate quintessence and phantom fields from non- viscous and viscous fluids in an anisotropic universe 

respectively. 

Quintessence and phantom fields are generally given by the action 



S 



d^x. 



1 



eiV^f ~ V{<t>) 



The energy density and pressure of scalar field (DE) are given by 



and 



4' + v{<t>) 



= -e(/.2-y(^). 



(36) 



(37) 



(38) 



where e = ±1. e = 1 is referred to as quintessence whereas e = —1 is referred to as phantom. From eqs. (|29|) . 
(|30|) and eqs. (|37|) . (|38|) we find the general form of the scalar field and potential V{(j)) as 



and 



==2e 



V{ct>) = 2 



Co, 



H\l + q) + La-^ - -H^n^'a-^ - ^VsQ^^ 



Co. 



ij2(l -q)- 4La~^ - -i?2o™a-3 + ^V3?]^ 



(39) 



(40) 



Note that putting ^o = and e = 1 in eqs. ([5^ . PO]) we get the scalar field and potential of quintessence. Also 
for sufficiently large time, the asymptotic behavior of (j) and V^cj)) is given by 



(-e^oVs)' 



t + constant. 



(41) 



and 



V{ct>) ^ CoVS, (42) 

respectively. Eq. (|^T|) clearly shows that the only possible scenario at far future is the phantom scenario as e = 1 
(quintessence) gives an imaginary (j). It is worth to mention that at late time i.e a — > oo which implies ^o ~^ 0, 
the potential asymptotically tends to vanish and (/) = constant. 



6 Statefinder Diagnostic 

V. Sahni and coworkers (2003) have recently introduced a pair of parameters {r, s} called ^' statefinders'^ , which 
are useful to distinguish different types of dark energy. The statefinders were introduced to characterize primarily 
fiat universe models with cold dark matter (dust) and dark energy. They were defined as 



r-n 



ii/3' 



3(-7 



On 
2 ' 



(43) 



Here the formalism of Sahni and coworkers is extended to permit curved universe models. If we suppose that dark 
energy does not interact with dark matter (as we assumed), then the statefinder pair can be further expressed 
as 

2 H 



= nr, 



9w'*^ 



(44) 



1 (') 

where il = 0„i + fJcje. The statefinder is a "geometrical" diagnostic in the sense that it depends upon the 
expansion factor and hence upon the metric describing space-time. 

If the dark energy is due to a scalar field the equation of state factor w'^'^ is given by 



02 + 2eV{(j)) 
by taking differentiation we get 



^^-2=?^ (^" )^-r +2^ r^-^i' (^1) 



where aij is the shear tensor which is given by 



Therefore, the statefinder s is also obtained as 



^' + 14 



3 



(46) 



^ 2e0(20y-02y) 
02 + 2ey(0) 
Using the equation of motion for the scalar field 

4> + 3H^ + eV = 0, (48) 

in eq. (|47| and inserting the result into (|44[) we obtain (note that V — V'4>, V — A ) 

-^ + ^^ + ^^ (49) 

Furthermore, from Raychaudhuri's equation 

- = ^Hip'T - ^/"^(l + 3^.'^^) - ip"(l + 3c.") - la^ (50) 

a 2 6 6 3 

we find 

^ _ ^0 „2T^l/or,de\T 2 2 , 1 1^)2 



cTij = Ui-j + -{ui-ku'^Uj + Uj-ku'^Ui) + -9{gij + UiUj). (52) 



(53) 



To study the behavior of viscous DE more precisely we consider a toy model in the next section. 

7 Test Model 

To examine our above general results we present a worked example in this section. For this propose we assume 
the following scale factor 

a{t) ^ sinh(i). (54) 

By assuming a time varying deceleration parameter one can generate such a scale factor (Amirhashchi et al. 
2011). It has also been shown that this scale factor is stable under metric perturbation (Chen and Kao 2001). 
In terms of redshift the above scale factor is 

^ z=^^-l. (55) 



1 + z ' sinh(i) 

In this case one can find the DE energy density p'^'^ , the bulk viscosity ^{p'^'^): deceleration parameter q, and 
average anisotropy parameter A^ as 

p-^^ = 3coth2(t) + 3Lsinh"^(t) - pj," sinh"3(t) 



= 3 



l + {l + zf 



+ iL{l + zf-p'^{l + zf 



(1 + Z)4 

,l + (l + z)2 



l-a 



3^1 



^l + {l + zf 

(1 + Z)2 



(1 + z)^ 
tanh^(t) 



+ 3L(l + z)''-p™(l + z) 



l + (l + z)2 



(56) 

(57) 
(58) 
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Figure 1 : The plot of the DE energy density p'^^ , average 
anisotropy parameter Am, and the bulk viscosity £,{p'^'^) 
vs. z for pg" = 0.24, L = 0.1, Co = 0.1. 



Figure 2: The plot of deceleration parameter q versus 
redshift (z). 
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hi 



^inh^ = 1(6? 
^^ + sinh2(t) 3^ 1 



&^ 



''\ + {l + zY 



(59) 



where AiJ^ = Hi — H{i — 1,2,3) and Hi — ^^ Hi = ^, and Hy, ~ ^ are the directional Hubble's parameters 
in the directions of x, j/ and z respectively. 

Figure 1 depicts the variation of energy density of p"**^. Am, and £,{p'^'^) versus redshift z. As it is expected 
all these parameters are decreasing functions and approaches to zero at late time (z = —1). The variation of 
deceleration parameter (DP) is also shown in Figure 2. From this figure we observe that the value of DP in 
present time is almost —0.6 which is in good agreement with the value of DP obtained from observations. Also 
at late time i.e z = —1, deceleration parameter tends to —1 as in the case of de-Sitter universe. 

By using eq. (jM)) in eqs. ([M)) . ([M]) . and (HOI) and after simplification the EoS of viscous dark energy wjj?, 
scalar field cj) and the potential V{4)) are obtained as 



w. 



1 / +2tanh^(i) + 1 + Lsinh"''(i) cosh"^(i) 
3 \\ + L sinh""* [t] cosh-2 {t) - il^ sinh"^ {t) 



yl — Oi 



Co- 



coth 



1-2q 



it) 



[Qde^a 



2e 



Co 



smh-'^{t) + Lfimh-^{t) - -r2™sinh"''^(i)cosh2(<) - — VSfT^ 



(60) 
(61) 



V 



2tanh-2(i) + l-ALsmh-^{t) - -Vf^ smh'^ {t) cosh^ {t) + — V3f7^ 



(62) 
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Figure 3: The plot of EoS parameter versus redshift (z) 
for f)^ = 0.24, L = 0.1. 



Figure 4: The plot of energy £7"* and fJ''^ versus redshift 
(z) for Vtl^ = 0.24, L = 0.1. 



We can re- write eqs. (pIH) - (|5^ in term of redshift as 



'■^'f 



1 2 I 7-_(i±£)' 

+ ^\+{\+zY ^ ^o"(i + ^y 



,1-a. [(l + z)-^+(l + z)- 



-2i1-2q 



V ^2 



= 2e 



1 



{vl'^'^Y 



^0, 



(1 -f zf + L(l + zf - ^r!^"(l + zf{\ + (1 + zf) - ^■^P^- 



1 



- 4L(1 + zf - \%\\ + zf{\ + (1 + zf) + ^Vsfl^ 



(63) 

(64) 
(65) 



(l + z)2 

The behavior of EoS parameter, lo'^^ , in terms of redshift z is shown in Fig. 3. It is observed that the EoS 
parameter is a decreasing function of z and the rapidity of its decrease depends on the value of ^o- We see 
that in absence of bulk viscosity the EoS always varying in quintessence region (red line/solid line) whereas in 
presence of viscosity EoS cross PDL and varying in phantom region. But at the later stage of evolution it tends 
to the same constant value i.e (jj'^'^ = — 1 independent of the value of ^o- This behavior clearly shows that the 
phantom phase i.e u)'^'^ < —1 is an unstable phase and there is a transition from phantom to the cosmological 
constant phase at late time. As we mention above, the phantom phase instability of the universe is because of 
the fact that the viscosity dies out as time is passing. 



The matter density fi™ and dark energy density ^'^^ can be easily calculated as 

f}" ^ f7'o"sinh"'''(t) = f7JJ'(l -F zf. 



n 



"^^ -l + L sinh-^(t) cosh-2(t) - 17™ sinh-3(t) = 1 + i-il±il-^ - f]™(l + zf. 



l + (l + z)2 



(66) 

(67) 



Also from above two equations we obtain the total energy density as 

fl-)-z)6 



n = ?l'" 



n 



de 



l + L 



l + il + zf 



The variation of density parameters Jl™ and ^'^'^ with redshift z have been shown in Fig. 4. Here, we observe 
that ft'^'^ increases as redshift decreases and approaches to 1 at late time whereas $7™ decreases as z decreases 
and approaches to zero at late time. 
For our model, the parameters {r, s} can be explicitly written in terms of cosmic time t or redsfift z as 
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Figure 5: The plot of il'^'^ versus 51™. The solid line in- 
dicates flat universe (L=0). The dots locate the current 
values of n"*^ and rj"* for L = 0, 0.05, 0.5 



s 0.3 
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Figure 6: s — r evolution diagram. The dots locate the 
current values of the statefinder pair {r, s}. 



tanh^(t) 



1 



l + il + zf 



and 



1 + L sinh"^(i) cosh"^(t) - tanh^(f) 

2 



[1 + 2 tanh2(t) + L sinh-''(t) cosh-^{t)] | [2 + L(l + z)2 (1 + (1 + z)^)] 



3 

2 



(69) 



(70) 



Figure 5 shows the values of Hq'^ and fi™ which are permitted by our model. From this figure we observe that 
for case L = which represents a spatially flat universe (57 = 1), J7q^ « 0.76 and 51™ « 0.24. These results are 
in good agreement with the CMB results, the supernova results, and the computed density of matter in clusters. 
Other models with L ^ 0, represent open universes with 57 < 1. 

Trajectories in s — r plane corresponding to different cosmological models are shown in figure 6. The dots in 
the diagram locate the current values of the statefinder pairs {s,r}. From this figure we see explicitly that the 
ingredient parameter L (or K) makes the model evolve along different trajectories on the s — r plane. It is 
worth to mention that the cold dark matter with a cosmological constant (ACDM) diagrams (spatially flat) 
corresponds to the fixed point {s, t}acdm = {0, 1}. From eqs. (|69|) and (|70|) we obviously see that {s, r} = {0, 1} 
at late time i.e z = —1. 



8 Concluding Remarks 

Phantom field models have been suggested in order to provide a theoretical support for the recent observation 
that mildly favor the EoS of DE crossing —1 near the past. A lot of studies have been done in this regard 
and many phantom field models have been proposed. Some of these models are evolving from quintessence to 
phantom called quintom. However, theses models suffer from two major problems i.e. (1) Instability of phantom 
field and (2) finite future singularity (big rip). In this paper we proposed a simple mechanism to alleviate 
these problems by introducing a special form of bulk viscosity i. e. H = —i^oH {p'^-'^Y in the cosmic fiuid. In 
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this mechanism first, viscosity causes dark energy which is varying in quintessence to pass phantom divided 
hne (PDL) and drop it to the phantom region but since viscosity is a decreasing function of time, as time is 
passing it dies out and w'^'' leaves phantom region and tends to — 1 at late time. Hence the problem of future 
singularity (big rip) does not occur in this scenario. To test the impact of the anisotropy parameter (L), we 
perform a statefinder diagnostic on this scenario. This diagnostic shows that the statefinder parameters can 
probe the anisotropy of the model. May be future SNAP would be capable of probing this effect. In summary. 
The general form of the EoS parameter of viscous and non- viscous dark energy has been investigated in this 
paper. It is found that the presence of bulk viscosity causes our universe to get to the darker region i.e phantom 
temporarily. It is worth to mention that since our anisotropic model behaves as isotropic FLRW universe at late 
time, as a result, the phantom does not survive in isotropic universe as well. Our results fulfil the theoretical 
requirement argued by Carroll et all (2003) which state that, to avoid the big rib problem, all phantom models 
should decay to cosmological constant at late time. Moreover, since we have not restricted our study to the 
maximally symmetric FLRW space-times, our results seems to be more general than those obtained on the bases 
of this isotropic universes. 
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